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Abstract 

We have measured the correlation function of Polyakov loops on the 
lattice in three dimensional SU(3) gauge theory near its finite temper- 
ature phase transition. Using a new and powerful application of finite 
size scaling, we furthermore extend the measurements of the critical 
couplings to considerably larger values of the lattice sizes, both in the 
temperature and space directions, than was investigated earlier in this 
theory. With the help of these measurements we perform a detailed 
finite size scaling analysis, showing that for the critical exponents of 
the two dimensional three state Potts model the mass and the sus- 
ceptibility fall on unique scaling curves. This strongly supports the 
expectation that the gauge theory is in the same universality class. 
The Nambu-Goto string model on the other hand predicts that the 
exponent v has the mean field value, which is quite different from the 
value in the abovementioned Potts model. Using our values of the 
critical couplings we also determine the continuum limit of the value 
of the critical temperature in terms of the square root of the zero tem- 
perature string tension. This value is very near to the prediction of the 
Nambu-Goto string model in spite of the different critical behaviour. 
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1 Introduction. 



Three dimensional SU(3) gauge theory has many properties in common with 
QCD. Lattice simulations of the theory show that the static quark potential 
is linear, implying confinement. They also give evidence for a mass gap and 
a nontrivial glueball spectrum [TJ. At finite temperature there is a phase 
transition to a state in which the energy density is approximately described 
by a gluon gas. It can therefore be expected that from this model one ob- 
tains important information about the mechanism of confinement and the 
deconfinement transition in QCD. Furthermore, there exist analytic approx- 
imations in three dimensional SU(N) gauge theory which predict the string 
tension and the glue ball spectrum at zero temperature [21 E]- These results 
are in quite good agreement with lattice calculations [HE]. It is therefore im- 
portant to measure the string tension at finite temperature and the value of 
the critical temperature, because the analytic calculations should eventually 
be extended to these observables. 

In a series of papers [HI El El E], we have investigated the SU(3) gauge the- 
ory at finite temperature in two spatial dimensions using lattice simulations. 
In [6] we have shown that in the high temperature phase above approximately 
1.5T C , where T c is the critical temperature, the theory can be dimensionally 
reduced to a gauge- Higgs model in two dimensions, by which one can give an 
excellent description of the long distance properties of the full theory. In [7] 
we have analysed the two dimensional gauge-Higgs model in great detail. In 
[8] we have investigated the thermodynamics of the three dimensional theory 
in the high temperature phase. We have shown, in particular, that the trace 
of the energy momentum tensor has a non-perturbative behaviour in a region 
above the phase transition, analogous to the results found in SU(3) in 3 + 1 
dimension [TO] and in full QCD HUfTJ. 

In [9] we studied the theory in the low temperature phase. By measuring 
the correlation function of Polyakov loops we obtained the finite tempera- 
ture string tension. We showed that it can be very well described by the 
Nambu-Goto string model, as was predicted in [TOJ [Ti] , but only up to a 
temperature T ~ 0.7T C , where T c is the critical temperature. One should 
remark, however, that there are analytic calculations for a general fluctuat- 
ing bosonic string, which in three dimensions give universal values for the 
terms in the expansion in (T/T c ) 2 up to (T/T c ) 6 [151 US E]- Not only do 
these terms, of course, coincide with those from the development of the for- 
mula in Nambu-Goto string model to this order, but their contribution is 
practically indistinguishable from the full model up to T ~ 0.75T C . However, 
such a short expansion cannot give any hint about a phase transition, while 
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the Nambu-Goto model predicts the existence of a critical temperature T c 
at which the string tension vanishes. Moreover, it gives a value for the non 
perturbative ratio T c j ^ l /ao, where <To is the zero temperature string tension. 
The ratio only depends on the number of transverse dimensions. It does not 
depend on the group SU(N). In fact, the Nambu-Goto string model gives a 
result for the approach of the finite temperature string tension to T c , which 
corresponds to the mean field exponent v = 1/2. From universality argu- 
ments it is, however, expected that the transition in the SU (3) gauge theory 
in two spatial dimensions is in the universality class of the two dimensional 
three state Potts model [H]. Support for this proposal has been found in 
lattice calculations [19J. 

Other comparisons of the Polyakov loop correlations in SU(3) in two 
spatial dimensions with the string model, in particular at lower temperatures, 
have been performed with a different technique in [20]. Recently, they have 
been extended to the group SU(6) [21]. 

In this article we report on an investigation of the model very near to the 
phase transition, a region which has not been studied up to now. First we 
calculate the critical couplings (3 C (N T ), where N T is the lattice extent in the 
temperature direction, extending earlier calculations [19, 23J to considerably 
higher statistics and larger lattices. For this purpose we introduce a new 
and powerful method based on the work by Binder on finite size scaling [22J . 
Further, using results from the literature for a^/oQ, where a is the lattice 
spacing we obtain from the continuum extrapolation a value for T c j y/ao, 
which can be directly compared with the prediction of the string model. 

We further investigate the correlation function of the Polyakov loops be- 
low but near the critical coupling. We extract the mass and the susceptibility 
for a large number of spatial lattice sizes and for a range of couplings near 
the critical one. We consider the finite size scaling functions of these observ- 
ables and find an impressive agreement with the universality class of the two 
dimensional three state Potts model. However, although the string and the 
gauge models thus have quite different critical behaviour, we find that their 
critical temperatures are extremely close to each other. 

The plan of the paper is as follows. The definitions relative to the 
Polyakov loop correlations, and the set up for a critical scaling analysis of 
the mass gap and the susceptibility are given in Section 2. In section 3, 
T c j g 2 is computed and the scaling properties of the susceptibility and mass 
gap as functions of T/T c exhibited. The ratio T c j \Zoq is finally computed 
and the gauge and string models compared. A summary and conclusions are 
proposed in a last section. 
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2 The Polyakov loop correlation function. Loop 
susceptibility and mass gap. 

A discussion of the simulations was given in [9]. We have used the same 
algorithm, which was ported to GPUs (Graphical Processing Units). In this 
section, we first give the formulae needed for our new analysis of the physical 
quantities of interest, and then recall the scaling properties expected from 
universality near the transition. 

2.1 The model and the correlation function. 

The action used in the simulations is the standard Wilson lattice gauge action 

S w (U,(x,r)) = - ^ReTrUp), (1) 

p 

where U^(x,t) denotes the group element on the link in the direction fi, 
whose origin is located at x = (xi, x 2 ) in space and at r in the temperature 
direction. The lattice has extension iVj x iV T . The variables (xi,X2,t) are 
defined on integer values, x\,X2 = 0,1, N$ — l and t = 0,1, N T — 1. The 
matrix Up is the product of link matrices around a plaquette. The constant 
(3 is the lattice coupling constant. Note that in this article we will never use 
(3 to denote the inverse temperature. The temperature and volume of the 
lattice are defined by 



f = aN T , (2) 
V = (aN s ) 2 . (3) 

where a is the lattice spacing. The coupling constant (3 in the action is related 
to the coupling constant g of the continuum action by 

hm(a/3) = - 2 . (4) 

In three space time dimensions, g 2 has dimension of energy, and can thus be 
used as the energy scale in the continuum theory. 

When Ns — > oo at fixed iVV there is a phase transitions with a critical 
coupling j3 c (N T ) to a state where the Z 3 symmetry of the theory is broken. 
This corresponds to a finite temperature transition into a deconfined state. 

To investigate the theory around the phase transition we performed nu- 
merical simulations in the neighbourhood of the transition for a large number 
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of values of the coupling constant (3 and lattice extensions N$ and N T . We 
study the order parameter, which in this case is the Polyakov loop. The local 
Polyakov loops are winding around the temperature direction: 

N T -l 

L(x)=Tr II U T (x,r). (5) 

r=0 

We define the projected correlation function G(z) between two loops at dis- 
tance z from each other in the x\ direction by 

G W = ^ r Y, Re (L(0)L*{x = (z,x 2 ))). (6) 

In our earlier work [9], devoted to couplings well below f3 c (N T ), we have 
represented this correlation function by 



G(z) = b cosh 



N s 

m z 

2 



(7) 



The subtraction of a disconnected part from G(z) would be required above 
(3 C (N T ). In ([2]), the coefficient b and the mass m were fitted to the data at 
z larger than some short distance cut off zq. This allowed us to extract the 
mass gap m and thereby the temperature dependent string tension 

^ = mN T . (8) 

Fits with Eq.dH) must give stable results with respect to changes of z Q , in 
which case they directly provide a reliable estimate of m, thus of the largest 
correlation length £ = 1/m. It corresponds to the simplest situation where 
the correlation in z decays as a pure exponential at large enough distances, 
and this was the case in the analysis which we reported in [9]. In such a 
case the correlation function can be as well represented in momentum space 
by a simple pole, in the statistical mechanics called the Ornstein-Zernike 
(OZ) behaviour. In general, the analytic behaviour in momentum space on 
the negative real axis of the square of the momentum is more complicated, 
involving, in the continuum limit, several poles and/or cuts. 

This is especially so if a continuous phase transition exists in the ther- 
modynamical limit Ns — > oo at (3 C (N T ). At the transition the correlation 
function is expected to decrease as a power in z, characterized by a critical 
exponent, the anomalous dimension 77. As a result, close to (3 C (N T ) and for 
values of (Ns, N T ) accessible in practice no simple parametrization of G(z) 
is available. Various parametrizations of its Fourier transform in the contin- 
uum limit have been proposed in the literature (for a review, see Ref. [25]). 
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They share the property that, apart from at /3 C (N T ), the nearest singularity 
remains an isolated pole in the complex plane of the momentum squared. 
Here, for any finite lattice, we will determine the mass gap squared as the 
distance to zero of the nearest pole. 

We start by transforming the numerical data for the correlation function 
to momentum space. Given G(z) on the integer values z — 0, 1, N s — 1, 



we define G q via 



where 



G(z) = J2e- iqz G g , (9) 



i 



27m , , 

g=-— ; n = 0,l,...,iV 5 -l. (10) 

l\ s 

Using the symmetries of the action, the periodic boundary conditions on the 
lattice and the reality of G(z), it is suitable to rename G q as G(p 2 ), where 
to q we associate 

P = 2sin(|). (11) 

Given numerical data for G(z) we compute the susceptibility x °f the 
Polyakov loop and the inverse of its Fourier transform H(p 2 ), 



Ns-l 



z=0 

X = JVfG(O), (13) 
H(P>) S ^ (14) 

For these measurements we restrict ourself to the disordered phase, where 
the subtraction of \(L)\ 2 is not necessary. On a finite lattice (L) = also in 
the ordered phase in simulations which are long enough because of tunneling 
between the degenerate vacua. 

Up to now, the functions G(p 2 ) and H(p 2 ) are defined only on the discrete 
values given in Eqs. fTTUl [TTjl . We now extend these functions to arbitrary 
complex values of p 2 , and define the mass m as the first zero of H(p 2 ) on the 
negative real axis, 

H(-m 2 )=0. (15) 
In the OZ approximation, we trivially get 

H z{p 2 ) = a + a ± p 2 , (16) 

m = — . (17) 
ai 
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Within our assumptions, the point p 2 = —m 2 , an isolated single pole of G(p 2 ), 
is generically (in the thermodynamical limit) inside the circle of convergence 
of the series expansion of H(p 2 ) around zero. Because on the one hand we 
focus on the long range properties of the correlation, and on the other hand 
wish to eliminate as much as possible discretization effects, we truncate the 
full series of H(p 2 ) to a small order n max in p 2 , and determine its coefficients 
from small p 2 data only. We write 

^max 

H{p 2 ) = *n( P 2 T, (18) 

n=0 

and solve for a n the system of n max + 1 equations 

an(p 2 iT = H{p 2 ) ; 2 = 0,l,---,n p (19) 

n=0 

where Pi are the discrete values defined by Eqs. ffTOj [TTj) and the right hand 
sides are measured using Eqs. (I11|I12|I14I) in n p points. Then m 2 is the smallest 
positive solution of Eq.([T5]). The case n max =l is the OZ approximation ffTol) ; 
the results presented in section 3 in the context of the scaling behaviour of 
m are those obtained with n max =2 and n p = 3. The errors quoted there 
are statistical errors only, estimated by a bootstrap technique applied to the 
whole set of configurations measured. 



2.2 The setup for the critical scaling analysis 

Near the phase transition one may expect that for given N T the theory is de- 
scribed by an effective Z3 symmetric two-dimensional model for the Polyakov 
loop ([5]), which is in the same universality class as the two dimensional three 
state Potts model [THj- It should therefore have the same two independent 
critical exponents v and rj as in the latter model. These are known from 
analytic calculations in that model to be 

v = 5/6, (20) 
r] = 4/15. (21) 

The universality hypothesis is supported by the results of [19]. There the 
critical exponents reported for N T = 2, 4 and 6 are, within errors, which are 
around 20% always compatible with the above expectations. 

Using the definitions (11211151) . we will assume finite size scaling behaviour 
for Ns large enough. Thus the mass gap and the susceptibility are represented 
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close to the transition by 




Ng m = 
X _ 



N^ s G(0) = f 2 (s), 



(22) 
(23) 



where 



s = N^(l 



W c (iv r )), 



(24) 



and 



fi(s -> oo) 
f 2 (s oo) 



'OC 



s 



s 



V 



u(2-rf) 



(25) 
(26) 



Here, (3 C (N T ), yet to be determined, denotes the critical value of the gauge 
coupling in the thermodynamic limit. The prefactors on the left hand sides 
of (I22II23I) are such that for (3 close to (3 C (N T ) the correlation length and 
the susceptibility behave in the N$ — > oo limit as prescribed by the critical 
exponents, that is 



In the following analysis, we assume universality to be true and use the 
values of v and r\ given in Eqs. fT20l |2~T|) . A posteriori the fact that the data 
fall on unique scaling curves show that these exponents are the correct ones 
for the SU(3) gauge theory in two spatial dimensions. For the susceptibility 
we make a further test of this assumption showing that the data do not fall 
on a universal function /^(s) if we assume the mean field exponent rj = 0. 

In the next section, we use finite size scaling applied to a different quan- 
tity, more easily measurable on the lattice, in order to determine (3 C (N T ) for 
N T = 4, 6, 8, 12. Then the scaling assumptions of Eqs. fT22l - 126]) are checked 
with a high accuracy in the case N T = 6, for which a large density of data 
points in parameter space has been acquired. Finally an extrapolation of 
Pc{Nt) to the continuum limit iV T — > oo is performed, which allows to com- 
pare the gauge theory and the string model near their transition points. 



m oc (l-P/p c (N T )Y, 
X oc (l-/3//3 c (N T ))-^\ 



(27) 
(28) 
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3 Applications of finite size scaling. Critical 
coupling, universal scaling functions and 
the continuum limit compared with the string 
model 

3.1 The critical couplings (3 C (N T ) 

We recall here the properties of finite size scaling, which we need to deter- 
mine the critical coupling. We determine /3 C (N T ) from the condition that the 
average value $^ of a classically dimensionless functional &{L} of the effec- 
tive field L(x) defined by Eq.flSD does possess finite size scaling properties 
analogous to those of m and x i n Eqs. (I22|l23p . Generically, we set 

N*9 L = f(s) (29) 

s = Ny v {i-p/p {N T )) 

f(s oo) oc s p , (30) 

where <ft, u, p are critical exponents. The conditions that, in a domain includ- 
ing = C (N T ), 

i) $l is well defined on the lattice for any Ns 

ii) $£, exists in the thermodynamic limit 
imply the relation 

p = v<p. (31) 

In this limit Eq. (j29~|) then gives 

$i oc (1 - 0/p c (N T ))"*, (32) 

and is the anomalous dimension of $. If (3 approaches (3 C (N T ), one may 
expand f(s) around zero and rewrite fl29|) as 

N* $i = a + bN 1 /^ - (3/f3 c (N T )) + (33) 

where a and b are constants. The right hand side is thus linear in N^ u at 
j3 fixed, and (3 C (N T ) is that value of (3 for which the slope 6(1 — (3/(3 c (N T )) 
vanishes. The method is, of course, particularly useful if <p and v are known. 

The order parameter (L(x)), which is classically dimensionless and has 
an anomalous dimension rj/2 in two dimensions, does not fulfill the condition 
i) above since, on the lattice, it vanishes even in the broken phase, due 
to fluctuations between degenerate groundstates. For the same reason the 
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susceptibility, which is quadratic in L and thus has the anomalous dimension 
rj (see Eqs. fl23ll26p ). is not suitable because its measurement via the definition 



(fi~3]) requires G(0) to be the connected correlation (by subtracting |(L)| 2 ). 
One could in principle use G(0) without subtraction, but we have found that 
it is less stable over the transition than the variable defined below. 
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Figure 1: Np /2 £ 3 

function of Ng for fixed values of ft. The temporal 
extent is iV r = 8. The straight lines correspond to the fits described in the 
text. 



In fact, we consider the quantity £3 defined as follows 

° x 

C 3 = Re(C 3 ). (35) 

This quantity is particularly interesting. C is the lattice average of L(x), and 
its cube is the monomial of lowest degree that does not suffer from the same 
disease as (L) itself: in the ordered phase £ 3 fluctuates around the same real 
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value whatever vacuum is chosen during the simulation. We have checked 
that the imaginary part of (£ 3 ) is always negligible, due to the reality of the 
action, and use its real part in ( |35|) for convenience. 
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Figure 2: The slopes of N s £3 for N T = 8 as a function of (3. The 
curve corresponds to the fit described in the text. 

In FigJT] we show data obtained for £ 3 at N T = 8 inside a domain in /3 
which strongly suggests that it contains a value of (3 such that £3 is indepen- 
dent of Ns- More precisely, inside the domain in {/?, Ns} shown, it appears 
that the data are consistent with £3 being linear in N§ , as expected if an 
expansion of the type of Eq. fl33l) is valid: The straight lines drawn actu- 
ally result from linear fits to the data at /3 fixed. They are restricted to go 
through the same point for vanishing argument, as demanded by Eq. (15^1) . 
Similar results are obtained for the other values of N T . In addition, the same 
ansatz (1331) implies that the slopes of £3 in N l J u at iVV fixed are proportional 
to 1 — /3//3 c (N T ) near (3 C . To take care of higher order corrections we fit the 
slopes, which we denote by £(/?), using 
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e(/3) = a 1 (/3-/3 c ) + a 2 (/3-/3 c ) 2 , (36) 

where a\ and a 2 are the parameters of the fit. These fits have very good 
X 2 /d.o.f. ~ 1. In FigfS] we show as an example the fit for N T = 8, which 
confirms the nearly linear behaviour of the slope near /3 C . The same is true for 
the other values of N T exploited. Those fits then complete the determination 
of /3 C (N T ) shown in the Table below. Note that the errors are the statistical 
ones only. The critical coupling /3 C (4) has been determined in two earlier 
works [19] [23]. The values were 14.74(5) and 14.717(17) respectively. The 
method in these investigations are different from each other and from ours. 
Systematic errors have not been taken into account. Therefore the results 
are in very good agreement with ours, giving further support to our method. 



N T 


C {N T ) 


4 


14.7404(7) 


6 


21.374(3) 


8 


27.952(4) 


12 


41.075(7) 



3.2 Finite size scaling for m and x 

Knowing f3 c (N T ), we are now able to test the universal scaling of the mass 
gap and the susceptibility. They are defined in Eqs.([22]-[26]). Their expected 
finite size scaling behaviours were described in Eqs.( TT5yT3l . Here we illustrate 
the fact that, apart from a known N$ dependent prefactor, they depend on 
the two parameters (3 and N$ via a function of the single variable s of Eq. 
(j24j) . We have measured the correlations, and thus m and \ f° r N T = 4 
and 6. Here we concentrate on N T = 6, our conclusions for N T = 4 being 
identical. We have used values for Ns from 48 to 264, and /3-values in the 
neighbourhood of /3 C = 21.374. 

We start with the mass because in this case only the critical exponent 
v enters the analysis and thus can be probed independently. The scaling 
properties of the mass described by Eqs. (12211241125]) are illustrated for N T = 6 
in Fig. [3j The quantity mNs is plotted versus s u /N T . The factor 1/N T is 
inserted in the latter variable because according to ([8]), the continuum limit 
of the product mN T depends only on the temperature, and thus provides the 
leading term (1 — T/T c ) v as iV T -> oo. The values of mN T collected in our 
earlier paper [9] for N T — 4 and 6 are, in fact, found to be close one to the 
other. 



12 



Up to s ~ 9 (s u /N T ~ 1), the existence of a finite size scaling law is 
impressively confirmed, showing a high density of points lying on one single 
curve. The curve is a fit to a second degree polynomial in s. This validates 
the value v = 5/6 expected from universality. For s > 9 the data in the figure 
lie in the neighborhood of a straight line, corresponding to the asymptotic 
behaviour of Eq. (I25p with the same value v = 5/6. 

We now probe the anomalous dimension 77 = 4/15 appearing in the left 
hand side of Eq.f l23|) . This is done in Fig. H]by plotting NgX against mNs- 
This choice eliminates the explicit dependence on the critical exponent u, 
and therefore is a direct test of rj. This successfully validates rj = 4/15. The 
curve in the figure corresponds to a fit 

^ (0 > = TTMF (37) 

where h = 2.53(2) and b 2 = 0.850(4). 

As already mentioned, to measure the mass gap and susceptibility above 
(3 C (s < 0) one should subtract from G(z) the disconnected contribution 
|(L)| 2 . On a finite lattice (L) vanishes, because of the tunneling between the 
degenarate vacua. In practise this happens in Monte Carlo simulations near 
the critical point for N$ not too large. A popular way out is to replace \ {L)\ 2 
by (\L\ 2 ). We will not use this procedure here. By the way, our choice to use 
£3 to identify the critical (3 was dictated by the need for a control parameter 
which evolves smoothly across the transition. 

The analysis described above, which has also been performed for N r — 4 
with the same conclusions convincingly comforts the expectation that the 
critical behaviour of the gauge model is governed by the same exponents as 
the two-dimensional three states Potts model. We have not tried to determine 
these exponents a priori from demanding a unique scaling curve for each of the 
variables considered, Cs,m,x- But, taking the susceptibility as an example, 
we show a contrario in Fig. [5] that choosing the mean field exponent 77 = 
is far from giving a unique scaling curve. 

3.3 The critical temperature in the continuum com- 
pared with the string model 

As shown above, for fixed N T there is a phase transition at a critical coupling 
fic{N T ) in the thermodynamic limit N$ —> 00. In the continuum theory there 
is thus a phase transition at the critical temperature T c , given from Eqs. (J21 
HDby 

-4 = lim ' c \ T Tj 38 
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N S /N T (1-T/T C ) V 



Figure 3: The scaling behaviour of the mass gap at N T =Q. The full line is 
the finite size scaling parametrization discussed in the text. The dashed line 
is the asymptotic behaviour given in Eqs. (|24f25p . 

To determine this ratio we use the data given in Table 1, and assume that 
the approach to the continuum limit is given by 

From the data for the largest values of N T , it is clear that the leading cor- 
rection is linear in 1/N T . Inserting the values from Table 1 in Eq. f !39j) . we 
obtain a quite good fit. We obtain 

^ = 0.5446(4). (40) 
9 

Another choice is to use the square root of the zero temperature string 
tension as energy scale. This is also what we need to compare our results 
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Figure 4: The scaling behaviour of the susceptibility at N T =Q. The curve 
corresponds to the fit described in the text. 

with the string model. In our earlier paper [9], we used data in the literature 
[23] [20] to obtain the following parametrization: 

a^^F{(5) = -^ (41) 

where 

h = 3.325(6) (42) 

z = 2.0(2) (43) 

p = 3.7(1) (44) 

The errors on the parameters are strongly correlated. What is important is 
the error on the interpolating function F ao {f5). We found that the absolute 
value of the error on l/F ao (j3) grows from 0.002 for — 14 to 0.01 for 
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Figure 5: This figure shows the absence of a unique scaling function for the 
assumption of mean field exponents. 

(3 = 42. The effect of this error on T c j ^Joq was estimated using the bootstrap 
technique. 

To obtain the continuum value of A /oo normalized to the scale setting 
parameter g 2 , we form the ratio 

= lim Bfl = 1 = o.554(l) (45) 

This value is in agreement with an earlier investigation which gives a value 
0.5528(8) for the same quantity [23], validating our parametrization of the 
data. 

The ratio T c j y/ao in the continuum limit is obtained from the ansatz 



N T F ao (f3 c (Nr)) v^ V iV r 2 N*' 
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We find 



0.986(1) 



(47) 



with x 2 /d.o.f = 0.07. This is lower than the value TJ ^ = 0.999(4) in Ref. 
[23] , where the same extrapolation formula is used, but for smaller values of 
N T = 2 - 5. 
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Figure 6: The quantity T c f \/oq as a function of 1/iVj?. The full line corre- 
sponds to the fit by Eq. (|46j) . The dashed line is the string model value in 
Eq. (ED}. 



As can be seen in Fig. 6 the data show clearly that the leading correction 
to the continuum value in f )46|) has to be quadratic in 1/N T as expected for 
a ratio between physical quantities. 

In the Nambu-Goto string model, the temperature dependent string ten- 
sion is predicted to be [14] 



a(T) = a \ 1 



y2 

c 



(4£ 
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where 




3(T 



(49) 



(D - 2)tt 



Here the constant D is the number of space time dimensions of the gauge 
model. The prediction only depends on the number of transverse dimensions 
and not on the group SU(N). In our case which has one transverse dimension 
it leads to 



Furthermore Eq. (JIB"]) shows that the approach to T c corresponds to the mean 
field critical exponent v — 1/2. 

It is quite astonishing that the measured value (pETj) for SU (3) is so close 
to the theoretical one fl50l) . in particular as we have shown that there is a 
scaling region, where the mean field behaviour of the string model is not 
valid. In principle it would be interesting to investigate if the prediction (150]) 
becomes better in SU(N) for larger values of N. It has, however, been shown 
that for N > 5 the transition is of the first order [2U1 |2"T] as is also predicted 
from the universality arguments [18]. Then a{T c ) has a jump, and Eq. ( 150]) 
is not expected to be valid. Values of T c /y/ao have also been obtained for 
SU (2) and Z 2 gauge theories. In these cases the values of this ratio is quite 
far from the string value in Eq. §U]), being 1.12(1) [27J and 1.222(2) [2B] 
respectively. 

4 Conclusions 

In this article we have presented results for the behaviour of three dimensional 
SU(3) gauge theory near the finite temperature phase transition. We have 
introduced a new and powerful application of finite size scaling to extract the 
value of the critical coupling. For this purpose, we have used the third power 
of the Polyakov loop, which evolves smoothly over the transition. To extract 
(3 C we assume the critical exponents of the two dimensional three state Potts 
model, which are expected to be valid because of universality arguments 
[18] . Employing this method, we have determined the critical coupling on 
considerably larger lattices both in space and temperature direction than in 
earlier investigations. 

That the phase transition is in the universality class of the abovemen- 
tioned Potts model is strongly supported by our finite size scaling analysis 
of the mass gap and the susceptibility. Their values were extracted from an 
analysis of the correlation function of the Polyakov loop in Fourier space. 
We used our measured critical couplings and these correlation functions at 




(50) 



18 



a large set of values of the spatial extension and the coupling constant near 
the transition. In this region all the data for the mass gap fall within errors 
on a unique scaling curve, if we assume the critical exponent v to be that 
of the Potts model. This test is sensitive to the exponent v only. A further 
check of the universality, which is sensitive to the exponent r\ was performed 
for the susceptibility in the disordered phase. Also in this case the data fall 
on a single scaling curve when we use the expected value i] = 4/15. This is 
not the case if we e.g. assume the mean field value 77 = 0. This analysis was 
performed for N T — 4 and 6 with consistent results. 

Finally we have used our values for the critical couplings for N T = 4, 6, 8 
and 12 and the zero temperature string tension determined at these values 
from an interpolation formula of data in the literature to make an extrapola- 
tion to the continuum and obtain the ratio of the critical temperature to the 
square root of the zero temperature string tension. We have compared this 
ratio to the corresponding result in the Nambu-Goto string model. Although 
the string model predicts a temperature dependent string tension which has 
mean field behaviour in the scaling region, the critical temperature in this 
model is, in fact, very close to the lattice result. The ratio is a non per- 
turbative and non universal observable. Thus it is a very important result, 
which may imply that the string model describes the spectrum and the mul- 
tiplicities of the excited states of the SU(3) gauge model in three space-time 
dimensions. 
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